Finite-Momentum Dimer Bound State in Spin-Orbit Coupled Fermi Gas 
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We investigate the two-body properties of a spin-1/2 Fermi gas subject to a spin-orbit coupling 
induced by laser fields. When attractive s-wave interaction between unlike spins is present, the 
system may form a dimer bound state. Surprisingly, under proper conditions, the bound state 
obtains finite center-of-mass momentum, whereas under the same condition but in the absence of the 
two-body interaction, the system has zero total momentum. This unusual result can be regarded as 
a consequence of the broken Galilean invariance by the spin-orbit coupling. Such a finite-momentum 
bound state will have profound effects on the many-body properties of the system. 
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With the experimental achievement of realizing syn- 
thetic gauge field in ultracold neutral atoms [1], there 
have been an intensive search for the emergence of ex- 
otic quantum phases [2-6]. So far, both Abelian and 
non-Abelian gauge fields have been realized for quan- 
tum gases using the two-photon Raman process [1, 7, 8]. 
The realized non-Abelian gauge field leads to a particu- 
lar spin-orbit coupling (SOC) which can be regarded as 
an equal-weight combination of Rashba and Dresselhaus 
SOC. Schemes of generating various other types of SOC 
have been proposed (e.g. [9-12]) and are under active 
experimental investigation. Previous theoretical studies 
have shown that the interplay between SOC and effective 
Zeeman field (induced by the same laser fields and/or real 
magnetic fields) can lead to intriguing quantum states 
with non-trivial topological properties. However the con- 
ditions under which such states can be realized have not 
been met in current experiments so far. 

In this Letter, we focus on the two-body properties 
of a degenerate spin-1/2 Fermi gas subject to spin-orbit 
coupling and effective Zeeman field. In the absence of 
two-body interaction, such a system exhibits zero total 
momentum. However, we show that in the presence of 
attractive s-wave interaction and with a proper combi- 
nation of the SOC and the Zeeman field, the system may 
form a dimer bound state with finite center-of-mass mo- 
mentum. At first sight, this result is very surprising 
since the s-wave interaction is manifestly momentum- 
conserving. A closer examination will reveal that this is 
a natural consequence of the broken Galilean invariance 
due to the presence of the SOC [13-15]. In principle, 
this finite-momentum dimer state can be realized in the 
current experimental system with equal- weight combina- 
tion of Rashba and Dresselhaus SOC. The momentum 
of the bound state can be measured through standard 
time-of-flight technique. 

We start by formulating the single-particle Hamilto- 



nian for a non-interacting homogeneous Fermi gas in 
three dimension: 
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where we have defined SOC strength vector v = 
(i>i,t>2,v3), and Zeeman field vector A = (Ai,A2,A3). 
v — (&xi <Jy,&z) are Pauli matrices acting on the atomic 
(pseudo-)spin degrees of freedom. This description is a 
general model valid for various coupling schemes. The 
single-particle spectrum can be straightforwardly ob- 



tained as — e k ± J X)Li( w «^i + ^i) 2 with e k 



h is. /(2m). The superscripts '±' define the two helic- 
ity basis which are related to the original spin basis by 
the transformation 
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Next we consider the attractive s-wave contact interac- 
tion between un-like spins which, in terms of the creation 
and annihilation operators for the original spin states, is 
represented by 



^ int C q/2+kt C q/2-ki C q/2-k'iC q/2+ k't ' 

kk'q 



(3) 



where g is the bare coupling strength to be renormal- 
ized using the s-wave scattering length a s . A general 
two-body wave function describing a dimer state with a 
definite center-mass momentum q can be written as [2[: 
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Inserting this wave function into the Schrodinger equa- 
tion ("Ho + 'Hint)|^ , )q = -Eq|'I , )q and after some lengthy 
algebra, we arrive at four coupled algebraic equations for 



the coefficients i/j acr i (k). Self-consistency requires that 
the energy of the dimer state, Eq, satisfies the following 
Cooper equation: 
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where £k.q = -Eq — e f +k — e f-ki and the interaction has 
been regularized as l/g = m/(4:Trh 2 a s ) — l/V l/(2ek). 
The explicit forms of ^ CT(T <(k) can also be found, which 
we list in the Supplemental Material [16]. 



We now apply this general formalism to the experimen- 
tally realized system whose single-particle Hamiltonian 
takes the following form: 
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Here ^ and ^J. are field operators for hyperfine spin 
states a. O is the two-photon Raman coupling strength 
and k r the recoil momentum of the Raman beams which 
is taken to be along the z-axis. Finally S represents 
the two-photon detuning. To get rid of the exponential 
terms, we introduce a local gauge transformation: 

^ t (r) = e- ifc -** t (r) , ^(r) = e ife "** 4 (r) . (6) 

Defining the spinor ?/> = (ip^, ipi) T , we can recast H as 



h 2 k 2 



2m 



k z cr z + —a x 
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where we have neglected a constant energy shift E r = 
h 2 k 2 /(2m) (the reoil energy) in Ho- One can imme- 
diately see that Ho above takes the form of Eq. (1) 
with v = (0, 0, h 2 k r /m) and A = (0/2, 0, S). 
The two-body interaction Hamiltonian takes the form of 
Eq. (3) after we define the momentum space operators via 
tp a — J2 k e lk ' r Ck a /VV, and is obviously invariant under 
the transformation (6). 



Before we discuss the two-body properties of the inter- 
acting system, let us briefly examine the non-interacting 
limit. The single-particle ground state occurs in the lower 
helicity branch, and is given by |p m i n — ) at momentum 
Pmin = (0, 0, hk ) with energy 
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where for notational simplicity we have defined h = 0/2 
and A = h 2 k r /m. The value of ko is obtained from 



aSmin _ H 2 k XjXkp + 5) 

dk ~ to y/h 2 + {\k + 5) 2 



0. 



Measured in the lab frame, the momentum of the particle 
prepared in this ground state can be calculated as 

Po = (Pmin + hk r z) sin 2 Pmin + (p min - hk r z) cos 2 6> Pmin , 

(10) 

where the two terms proportional to hk r are included 
to 'undo' the local gauge transformation performed in 
Eq. (6). After some straightforward algebra, we can show 
that Po = 0, i.e., particle in the ground state has exactly 
zero momentum in the lab frame. 
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Figure 1: (Color Online) (a) and (b) Single-particle dispersion 
in the lower helicity branch. The horizontal dashed lines rep- 
resent the chemical potential of the Fermi sea. (c) and (d) Mo- 
mentum distribution along the z-axis for the non-interacting 
Fermi gas represented by (a) and (b), respectively. The rela- 
tive population of the two spin components are listed in the 
plot, blue dashed line for spin-down atoms and red thick line 
for spin-up atoms, (e) and (f) Momentum distribution along 
the z-axis of the dimer bound state, (g) and (h) Difference 
between n-|-(fc z ) and k z ) for the cases shown in (e) and 
(f), respectively. For the left column [(a), (c), (e), and (g)] 
we have 5 = and h = 0.6£V; for the right column [(b), (d), 
(f), and (h)] we have S = 0AE r and h = 0.6E r . For the 
interacting case [(e)-(h)], the interaction strength is given by 
l/(k r a 3 ) = 1. 



bound state only occurs on the BEC side of Feshbach 
resonance with a s > [17]. More importantly, qo = qqz 
will be non-zero and along the z-axis as long as both h 
and 8 are finite. Figure 2(b) displays go as functions of 
h and S. We note that qo is an even function of h, and 
an odd function of 5. 




Figure 2: (Color Online) Binding energy et (in units of E r ) 

(a) and the magnitude of the momentum qo (in units of k r ) 

(b) of the dimer bound state, as functions of Zeeman field 
strengths S and h. The white region is where no bound states 
can be found. The scattering length is given by l/(k r a 3 ) = 1. 



Next we consider a filled Fermi sea at zero tempera- 
ture. Instead of taking p = p m i n as in Eq. (10), we need 
to sum over all the momentum under the Fermi surface. 
Here again one can show analytically that the total mo- 
mentum of the Fermi gas is zero in the lab frame (see 
Supplemental Material [16]). The single-particle energy 
dispersion and momentum distributions of two examples 
of non-interacting Fermi sea (one for S — and the other 
5 7^ 0) are illustrated in Fig. l(a)-(d). 

Now we are in the position to discuss the two-body 
results when interaction is included. Following the gen- 
eral formalism, for a given set of parameters h, 6 and a s , 
we can obtain numerically the eigen-energy of the dimer 
state Eq as a function of q. The momentum qo at which 
Eq reaches the minimum labels the ground dimer state. 
The binding energy is defined as 

£{, = 2E min — Eq . (11) 

Only when ej, > 0, can we consider the dimer as a truly 
two-body bound state. Otherwise, its energy lies in the 
single particle continuum. Figure 2(a) shows £(, decreases 
with increasing h and S. Beyond a critical boundary 
value, binding Energy becomes negative and no stable 
bound state could be found. For this system, two-body 



For the two-body wave function given in Eq. (4), the 
momentum distribution of the hyperfine spin states in 
the lab frame is given by [16] : 

»t(k + M) = |^n(k-q/2)| 2 + |^tt(k-q/2)| 2 (,12) 
*H(k-M) - l^t(k-q/2)| 2 + |^ U (k-q/2)| 2 (.13) 

Both of these distribution functions will be symmetric 
along the x- and y-axis, which yields the average momen- 
tum along these two axis P* = J dk k x ^ y n a {k) — 0. By 
contrast, n a (k) will be in general asymmetric along the 
z-axis which results in finite values of = J d\£k z n a (k). 
The total momentum can be shown as [16] 

Plab - Pf ab + Pf ab = q + (JV t - NJkrZ , (14) 

where N a — J dkn cr (k) is the polulation in spin-er and 
satisfy the obvious constraint iV-f + N± = 2. In Fig. 1 (e) 
and (f), we illustrate how n a (k z ) — J dk x J dk y n a (k) 
changes without and with detuning S. To see it 
more clearly, we plot the difference between n^{k z ) and 
n^(—k z ). Notice that it is SO coupling that breaks spatial 
reflectional symmetry such that n a {k) ^ n a {— k) with 
er =t-|- However, for <5 = 0, one still has the symme- 
try %(k) = n±(— k) as measured in the experiment of 
Ref. [7]; with finte 5, this symmetry between n-f(k) and 
n_i(— k) is further broken. 
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Note that for either h — or S — 0, we obtain q = 
and Nf = N± — 1. In either of these case, both P\ and 
P£ will be finite, but they have equal magnitude and 
opposite sign and hence the total momentum P z = 0. 
When both h and S are non-zero, we obtain finite qo and 
moreover iV-j- =/= JVj.. Our numerical calculation shows 
that, under such circumstances, P 2 w 2qo- 




Figure 3: (Color Online) (a) For spherical SOC, bound state 
energy E q as a function of q (momentum along the z-axis) for 
different values of Zeeman field strength, from top to bottom 
5/E r = 0,0.2,0.4,0.6,0.8 and l/(k r a s ) = 1/2 for all curves. 

(b) Ground state momentum qo as a function of 8 at different 
values of scattering length. From left to right, the curves cor- 
respond to l/(k r a s ) = —I, —1/2, 0, 1/2, and 1, respectively. 

(c) go as a function of l/(k r a s ) at S/E r = 0.1. 

This situation, combined with the results for the non- 
interacting system, leads to the following peculiar phe- 
nomenon: In the absence of two-body interaction, the 
system possess zero total momentum; when the inter- 
action is turned on, the system picks up a finite total 
momentum even though the interaction Hamiltonian (3) 
seems to be momentum-conserving. This peculiarity is 
indeed a salient feature of the SOC. Since the linear mo- 
mentum of the atom is intimately coupled to its (pseudo-) 
spin, as the interaction induces redistribution of atomic 
population (i.e., change the value of a z ) as shown in 
Fig. 1(d) and (f), it also modifies the total momentum of 
the system. This phenomenon can be regarded as a man- 
ifestation of the broken Galilean invariance, which is also 
responsible for other unusual behaviors such as the de- 
viation of dipole oscillation frequency in a harmonically 
trapped system [13, 14], and the ambiguity in defining 
Landau critical velocity in spin-orbit coupled condensates 
[15]. 

The momentum of each spin species can be readily 
measured in experiment using the time-of-flight tech- 
nique. This has actually been performed in a non- 
interacting SOC Fermi gas [7]. However, our calculation 
shows that qo is only on the order 1% of recoil momentum 
for the equal Rashba-Dresselhaus SOC (see Fig. 2) which 
is the only type of SOC realized so far. The magnitude 
of qo, however, can be greatly enhanced under other SOC 



scheme. As an example, we consider the spherical SOC 
coupling scheme that is recently proposed [11, 12]. For 
this case, the SOC strength v\ = t>2 = «3 = v. We take 
the spin-orbit coupling strength v — h 2 k r /m to be the 
same as in the previous case. Due to the isotropic na- 
ture of the SOC term, the direction of the Zeeman field 
is irrelevant. We choose it to be along the z-axis, i.e., 
A = (0, 0, S), which, as we shall show, leads to a dimer 
state with finite moment qoz along the z-axis. Fig. 3(a) 
demonstrates that, with increasing S, the minimum of 
of the bound state energy E q deviates away from zero to 
some finite value. Compared to the former equal Rashba- 
Dresselhaus SOC case, an essential difference here is that 
a two-body bound state can be realized even when a s < 
[17]. Fig. 3(b) and (c) show how q varies as functions of 5 
and a s . As one can see, as long as 8 is non-zero, the dimer 
ground state possesses finite momentum. Furthermore, 
The magnitude of qo can reach as high as 0.2k r . Such a 
large value should be easily detected in experiment. 

Finally, for the sake of completeness, we comment 
on the case of Rashba SOC with v = (vi, t>2, 0) and 
vi = V2 = h 2 k r /m. In this case, when an in-plane Zee- 
man field (i.e., with a component in the x-y plane) is 
present, the resulting dimer bound state will again have 
finite momentum. A plot similar to Fig. 3 can be ob- 
tained. Our calculation shows that the maximum qo can 
be reached is about 0.05fc r . 

In summary, we have studied the two-body properties 
of a degenerate spin-1/2 Fermi gas subjected to spin-orbit 
coupling and effective Zeeman field. We show that, two 
fcrmions, via s-wave scattering, may form a dimer bound 
state with finite center-of-mass momentum, under proper 
configuration of SO coupling and Zeeman field. We at- 
tribute this peculiar phenomenon to the broken Galilean 
invariance induced by the SO coupling and special role 
played by the Zeeman field. We have directly solved the 
two-body Schrodinger Equation in the general formalism, 
and considered momentum distribution in the lab frame 
for the experimentally realized system, where finite but 
relatively small bound state momentum go is found. Fi- 
nally, we demonstrate that the recently proposed system 
with spherical SOC [11, 12] can result in dimer bound 
state with up to 20% of the recoil momentum. In a many- 
body setting, our two-body calculation is relevant in the 
limit where most of the atoms form tightly bound pairs. 
When this limit is not reached, we still expect that the 
interplay between SOC and Zeeman field may lead to ex- 
otic superfluid state. For instance, as a direct analog of 
the finite-momentum dimer, the many-body system may 
support finite-momentum Cooper pairs [18] reminiscent 
of the Fulde-Ferrell state studied in the context of Fermi 
gases with spin imbalance. A systematic study of such 
many-body properties will be conducted in the future. 
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In this supplemental material, we present detailed expressions of the two-body wave-function and discuss about 
momentum distribution for both the non-interacting system and two-body bound state. 



A. GENERALIZED COOPER EQUATION 



First, we formulate the single particle Hamiltonian Eq. (1) (main text) in terms of creation and and annihilation 
operators, and denote S(k) = + A3 + v 3 k z and T(k) = — A 3 — v 3 k z , U(k) — v\k x + Ai — i(v 2 k y + A 2 ) , 
V(k) = v\k x + Ai + i(v 2 k y + A 2 ), then 



n = J2 (s( k ) 



(15) 



The lengthy part of the derivation of Eq. (5) (main text) concerns about simplifying - Ho|*) q . We proceed as the 
following, 

H i*) q = $>t4(k)[s(| +k)4 +kt +v(|+k)4 +k j4_ ki 
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+^ t (k)[s(| - k )4_ kt + v q k) C |_ k j4 +H (-i) 

+^rt(k)[S(| + k )4 +kt + + k)4 +k Jc|_ kt 

+^n(k)[s(3 - k )4_ kt + i/(f - k)4_ k j4 +kt (-i) 
+^(k)[T(| - k )4_ H + ej(S - k)4_ kt ]4 +kt (-i) 
+^ t (k)[T(5 + k )4 +H + [/(| + k)4 +kt ]4_ kt 

+Vu(k)[T(| + k)4 +H + [/(| + k)4 +kt ]4_ ki 

+Vu(k)[T(| - k)4_ H + Uq k)4_ kt ] c t f+H (-l) 

= E' + k)^n + - k)Vtt + r (§ - k )V>n + ^ (§ + k)^u 
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+ - k)^ t + n§ + k)Vn + r(f + k)^t + c (§ - k)^ u 



+ 



S(3 + k)^ t + (§ - k)^n + S(5 - k)^t + l/(| + k)^t" 

T(5 + k)V> U + - k)^ t + r(f - k)^u + + k)Vn' 
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"■("k"!" ^ kj, 
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"2* — I - k,J^ kj^ 



where J^ k denotes summation over positive momentum fc z . While interaction term is relatively simpler, 



Wint|*>q=fE'K( k )-^t(k) 



f t 
2"~| _ k'|~ 2" kj. 



After some simplifications, we obtain four coupled equations as 

(£k, q - 2« 3 fcz)^(k) = fE'KlkO-^tlkOj + fnf -k)]v n (k)+[t/(| + k)]^ u (k) 

k' 

fE'K(k')-^u(k' 



(£k,q + 2w 3 fc z )V'4,t(k) 
(£ k ,q - 2Ai - u 3 fe z )^ tt (k) 

(£k,q + 2A 3 + W 3 fe)Vu(k) 



y(| + k)]^ t (k) + [c/(|-k)]v^ u (k) 



[[/(|-k)]^ Ti (k) + [[/(| + k)]^ t (k) 
+ k)]vti(k) + [y(|-k)]^ t (k) 



where £ k 



q - £q 



Let us introduce 



3 the singlet wave function V's(k) = ^[V'T4.(k) — V\|.t(k)] f 



(16) 

(17) 

(18) 
(19) 
and the 



triplet wavefunctions ip t (k) = ^[V'u(k) + ^t( k )], ipu = 75 (V'Tt(k) + ^u( k )), ^ = 75 (V'tt(k) - ^u( k )), and after 
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adding and subtracting equations, we have 

£k,qipt ~ 2v3k z 4> s = (viq x + 2A 1 )tl> u + i(v 2 q y + 2A 2 )i> v 
^k.qV's - 2v 3 k z ip t 



2<? 



V 



^ Vs(k') - 2vik x tp v - 2v 2 ik y ip u 



£k,qi>u - (2A 3 + v 3 q z )ip v 
£k,qipv - (2A 3 + v 3 q z )ip u 



(v\q x + 2A 1 )tp t + 2iv 2 k y ip s 
-2vik x ip s - i{v 2 q y + 2A 2 )ip t 



(20) 

(21) 

(22) 
(23) 



or, in matrix form, 



( -2v 3 k z £ k>q -(v 1 q x + 2A 1 ) -i(v 2 q y + 2A 2 ) 

fk, q -2w 3 fc z 2w 2 fc y 2v 1 k x 

-2iv 2 k y -(viq x + 2Ai) £ kiq -(2A 3 + w 3 g' z ) 

\ 2«ife a: i(w 2 % + 2A 2 ) -(2A 3 + v 3 q z ) £ k)q 



As a first step, we focus on the subsequent sub-matrix equation 

£k, q -(v 1 q x + 2A 1 ) ~i(v 2 q y + 2A 2 ) 

-(v 1 q x + 2A 1 ) £ k , q -(2A 3 + v 3 q z ) 

i(v 2 q y + 2A 2 ) -(2A 3 + v 3 q z ) £ kjq 



\ 


<*,\ 




( ° \ 








f EkV- s 











J 








it' 




f 2v 3 k z ip s \ 


tpu 


H 


2iv 2 kyip s 


i>v , 




v -2vik x i> s J 



Denote 

M E 



£k, q -(wife + 2Ai) -i(v 2 q y + 2A 2 ) 

~{v 1 q x +2A 1 ) £ k , q -(u 3 (7 z + 2A 3 ) 

«(^2% + 2A 2 ) -{v 3 q z + 2A 3 ) f kiq 
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Using Cramer's rule in linear algebra, we express triplet components in terms of V's(k) as 
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<Mk) = 
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2^ s (k) 
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2^ s (k) 
M 



w 3 fc 2 -(uife + 2Ai) -i{v 2 q v + 2A 2 ) 
iv 2 k y £ ki q -(w 3 g z +2A 3 ) 

-uifea; -(v 3 <? z + 2A 3 ) £ k , q 

£k, q «3fc z -i{v 2 q y + 2A 2 ) 

-(wife + 2Ai) w 2 fc y -(w 3 q z +2A 3 ) 
«(«2% + 2A 2 ) -v\k x £ k , q 

£k, q -(«ife + 2Ai) v 3 k z 
-(«ife + 2Ai) £ k>q w 2 fc y 
i(v 2 g y + 2A 2 ) -(w 3 g z + 2A 3 ) -uifcj. 



(24) 
(25) 
(26) 



Then we insert expressions of V't(k), V'u(k) and ^>„(k) into Eq. (21), and integrate both sides over momentum and 
divide by the constant y Ek V's(k), we have the generalized Cooper Equation, 
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£ k , q (£^ q -ELiK'J l +2A,) 2 ) 



(27) 



and the singlet wavefunction is determined from Eq. (21) accordingly. The unnormalized singlet wavefunction is given 
by 



V'a(k) 



^k,q — 



4^, q E-=i vM - 4 (Eli vMvtqt + 2Ai) 



^,q (^, q -Ei=l(^ + 2A i ) 2 ) 

and the triplet wavefunctions are obtained from Eqs. (24), (25), and (26), correspondingly. 



(28) 
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B. SINGLE PARTICLE MOMENTUM OF THE GROUND STATE IN LAB FRAME 

For a single particle described by Hamiltonian (8) in the main text, in order to calculate the center-of-mass mo- 
mentum in lab frame, we need to "undo" the local gauge transformation, from spinor ^) to $. First, we notice that 
the ground state minimum stays in the lower helicity branch, which is a superposition of up and down components, 
|k— ) = — sin#k|k t) + cos#k|k I) and the minimum occurs at k m j n = (0, 0, fc ) where (h = m = 1) 

A(Afc + -5) 



While, from bare spin basis ( atomic pseudospin basis ) to the model Hamiltonian's basis, spinors are related by 



*,(r) = e-**^(r) 
For the ground state with quasimomentum k m j n , the momentum in the lab frame is 
P Z = |sin6» kmin | 2 (fco + fc r ) + |cos0 kmi J 2 (fco-fcr) 



(29) 

(30) 
(31) 

(32) 



1 L Xk x + S 



, \ 1 ( Xk x + 5 
(k + K) + - 1 + 



(k - k r ) 



(Afc + 5)k r k r 
= fc ; = k a rkn = 

where in the last line, Eq. (29) and k r = A in the dimensionless unit are used. 

C. TOTAL MOMENTUM OF NONINTERACTING FERMI GAS AT ZERO TEMPERATURE 

Now consider a non-interacting Fermi sea described by Hamiltonian (8) in the main text. For a given chemical 
potential /z, Fermi surface is defined by the following equation 

k ^ ~\~ k | 



^K 2 + (Xk z + S) 2 = /i 



(33) 



where we have assumed that fi lies below the bottom of the upper helicity branch and the Fermi surface is simply 
connected, as illustrated in Fig. 1(a) and (b) of the main text. The total mechanical momentum in lab frame is 
obtained by integrating over the three-dimensional volume of Fermi sea: 



Kb = Yl [\ 9 k \ 2 (k z + k r ) + \ cos 6 k \ 2 (k z -k r )] 
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dk 7 



k z 



(Xk z + 5)k r 
^h 2 + (Xk z + S) 2 



fi - -k 2 z + y/h? + (Afc, + S) 2 



V 



8(27t) 3 




(hi -kf + 8SX(ki - k 2 ) - 4(/i + \ 2 ){k 2 2 - k\)) 



where k\^ are the intersections of the Fermi level with the lower helicity branch at the z-axis and are simply the roots 
of Eq. (33) after taking k± = 0. Similarly, 



^lab = J2 [I Sin ^| 2 fcx + | COs6 k \ 2 k x ] = 
k 

P n = [\ sin e k \ 2 k y + \ cos e k \ 2 k y ] =0 
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D. MOMENTUM DISTRIBUTION OF THE DIMER BOUND STATE IN LAB FRAME. 

After identifying the ground state dimer momentum q , at which E q reaches the minimum, we consider the mo- 
mentum measured in lab frame. First, the density of spin-up atom at momentum p + k r z is given by 



n t (p + M) = (*|c^c pt |*) 



(34) 



iMp --*■)- M-z- - p) c n„-pi + [Mp - tt) - iM-r- - p) c , 



,qo 



qo-pt 



10) 



(2^ n (p - |)< o _ p , + a^tCp - f )<o-pt) 1°) 



= IV>n(p-f)l 2 + IV<TT(p- f)l 2 



and similarly the density of spin-down atom at momentum p — k r z is given by 

2 



n+(p - M) = (*|4c p J*) = |^ t (p - ^)| 2 + |V U (P - f )| 2 



(35) 



(36) 



The total momentum as measured in the lab frame is therefore 

Plab = J2 IP( n t(P) + %(P))] = qo + (N t - A^)M 



where N a = ?v(p) is the population of the spin-cr component. 



